Astronomy <fe Astrophysics manuscript no. H4065 


February 2, 2008 


(DOI: will be inserted by hand later) 





Spectral and Photometric Evolution of Young Stellar 
Populations: the Impact of Gaseous Emission at Various 

Metallicities* 



Peter Anders, Uta Fritze 



Alvensleben 



m 

O 
O 

(N 

X> 



> 

(N 
O 

m 
o 

a,: 

6 

CO 



Universitats-Sternwarte, 
ufritze@uni-sw.gwdg.de 



Geismarlandstr. 11, 37083 Gottingen, Germany; Email: panders@uni-sw.gwdg.de, 



X 



Received xxx / Accepted xxx 

Abstract. We include gaseous continuum and line emission into our GALEV models for the spectral and photometric 
evolution of Simple Stellar Populations (SSPs) for various metallicities in the range 0.02 < Z/Zq < 2.5. This 
allows to extend them to significantly younger ages than before. They now cover the age range from 4 Myr all 
through 14 Gyr. We point out the very important contributions of gaseous emission to broad band fluxes and 
their strong metallicity dependence during very early evolutionary stages of star clusters, galaxies or subgalactic 
fragments with vigorous ongoing star formation. Emission-line contributions are commonly seen in these actively 
star-forming regions. Models without gaseous emission cannot explain their observed colors at all, or lead to wrong 
age estimates. We use up-to-date Lyman continuum (= Lyc) emission rates and decided to use recent empirical 
determinations of emission line ratios relative to H/3 for subsolar metallicities. We justify this approach for all 
situations where no or not enough spectral information is available to determine all the parameters required by 
photoionization models. The effects of gaseous line and continuum emission on broad band fluxes are shown for 
different metallicities and as a function of age. In addition to the many filter systems already included in our 
earlier models, we here also include the HST NICMOS and Advanced Camera for Surveys (= ACS) filter systems. 

Key words, globular clusters: general, open clusters and associations: general, Galaxies: star clusters, Galaxies: 
evolution 



1. Introduction 

In Schulz et al. (2002) we presented spectral evolution- 
ary synthesis models for SSPs of a wide range of metal- 
licities 0.02 < Z/Zq < 2.5. Our models are based on 
isochrones from the Padova group that include the TP- 
AGB phase and on the spectral model atmosphere li- 
braries from Lejeune et al. (1997, 1998), rectified to yield 
agreement with observed colors from U through K of stars 
with effective temperatures in the range 2 000 K through 
50 000 K. Stars with higher temperatures are treated as 
black body radiators. These models give the time evolu- 
tion of spectra as well as of luminosities and colors in a 
large set of filter systems. As they did not include gaseous 
emission important during the lifetime of massive ionising 
stars, they cover a range of ages from 140 Myr through 14 
Gyr. 

SSP model results are useful not only for analyses of 
star clusters, the genuine SSPs, but are also readily su- 
perposed to describe the evolution of composite stellar 



* Data files are provided in the electronic version/at CDS, 
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populations like galaxies with star formation histories ex- 
tended in time and with various chemical enrichment his- 
tories. They can also directly be combined with cosmo- 
logical structure formation scenarios that include a star 
formation criterium. 

Systems with active ongoing star formation like very 
young star clusters, actively star-forming or star-bursting 
galaxies in the local as well as in the high-redshift 
Universe, protogalaxies and subgalactic fragments in par- 
ticular, require the inclusion of gaseous emission. This is 
not only important for the characteristic emission lines in 
the spectra but gaseous emission, both in terms of lines 
and continuum, give important contributions to broad 
band luminosities and colors (cf. Kriiger et al. 1995, 
Zackrisson et al. 2001). 

The wealth of HST ACS data currently becoming 
available - on many very young star cluster systems in 
particular - prompted us to not only include the gaseous 
emission in an updated way into our spectral evolution- 
ary synthesis models for SSPs of various metallicities but 
also to provide versatile tables for the luminosity and color 
evolution of these SSPs in ACS filter bands. 
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Results of these models have already been proven very 
useful in the interpretation of the ACS Early Release 
Observations of the very young star cluster systems in 
the Mice and Tadpole galaxies (de Grijs et al. 2002). 

2. Input Physics 

2.1. SSP Models 

Except for the gaseous emission that was not yet in- 
cluded in our SSP models before, we use the same 
input physics as presented in detail in Schulz et al. 
(2002). This includes isochrones from the Padova group 
containing the TP-AGB phase and model atmosphere 
spectra from Lejeune et al. (1997, 1998), extending 
from 90 A through 160 /im, for five different metal- 
licities Z = 0.0004, 0.004, 0.008, 0.02 = Z Q and 0.05 or 
[Fe/H] = -1.7, -0.7, -0.4, and +0.4. Inclusion of 
the TP-AGB phase has been shown to be very impor- 
tant for colors like V — I and V — K. As we have shown 
in Schulz et al. (2002), age-dating of star clusters on the 
basis of their V — I colors, as often done for young cluster 
systems in interacting galaxies or merger remnants, can go 
wrong by a factor 2 (at fixed metallicity and extinction) 
if the TP-AGB phase is neglected in the models. 

We present our results for a Salpeter and a Scalo IMF 
as in Schulz et al. (2002). 

2.2. Gaseous Emission 

Gaseous emission is primarily related to very hot stars, i.e. 
to massive stars in early evolutionary phases of an SSP. 
However, hot white dwarfs (WDs) can also contribute in 
later stages (but see next subsection). Gaseous emission 
depends on metallicity in a two-fold way. First of all, stars 
get brighter and hotter on average in stellar populations 
at lower metallicity. The lifetimes of low mass stars get 
shorter, those of high mass stars get longer at lower metal- 
licitics as compared to solar. This affects the output rates 
in terms of hydrogen ionising photons, Nj_, yc , of a stellar 
population as well as its absorption line spectrum and, 
hence both the strengths of all gaseous emission lines in a 
uniform way and the gas continuum emission. Second, the 
chemical composition and the physical properties of the 
gas exposed to the same ionising radiation field determine 
the relative strengths of different emission lines. 

The physical properties and the chemical composition 
of the gas ionised by a bunch of stars of given metallicity 
are not known a priori. They determine, however, the flux 
ratios of non- hydrogen element lines relative to H3. We 
therefore chose to assume that the gas has the same solar- 
scaled abundances as our single burst single metallicity 
stellar population and to use observationally determined 
emission line ratios for all non-hydrogen lines. In this re- 
spect our approach is similar to that of Perez-Gonzalez et 
al. (2002) who also use observed line ratios for their sample 
of strong H Q emitting UCM survey galaxies. It differs from 
those of Chariot & Longhetti (2001), Moy et al. (2001), 



or Zackrisson et al. (2001), who assume electron densities 
and temperatures for the gas at low metallicities and cou- 
ple a photoionisation code to their evolutionary synthesis 
model. 

On the basis of effective temperature and bolometric 
luminosity (and hence radius and surface gravity) of ev- 
ery star in a given isochrone, its flux of hydrogen ionising 
photons (A^Lyc) is calculated from up-to-date non-LTE ex- 
panding model atmospheres that take into account line- 
blanketing as well as stellar winds and recent tempera- 
ture and gravity calibrations (Schaerer & de Koter 1997, 
Vacca et al. 1996, Smith et al. 2002). Summing up the 
A^Lyc of all stars present in one isochrone gives the total 
A^Lyc of the stellar population at a given time. Depending 
on the metallicity a fraction of this flux is absorbed by 
dust immediately (30% is assumed to be absorbed for 
metallicities > 0.008, no absorption is applied for lower- 
metallicity environments, following Mezger 1978, see also 
Weilbacher et al. 2000). From the remaining A^yc flux of 
an isochrone we calculate the gaseous continuum emission 
and the hydrogen line fluxes as described in Kriiger et al. 
(1995) and Weilbacher et al. (2000) assuming ionisation- 
recombination equilibrium (Osterbrock case B). The de- 
tailed formulae were already presented in Kriiger et al. 
(1995). For instance, the A^yc is converted into an H3 
flux using 

F(U P ) = 4.757 x HT 13 • N Lyc . 

Emission line fluxes for elements other than H are calcu- 
lated from the line ratios relative to H^, given in Table 
1. 

For low metallicity gas, line ratios for an exhaustive 
set of strong forbidden and allowed transitions from UV 
through NIR are obtained from the extensive observa- 
tional database of Izotov et al. (1994, 1997) and Izotov & 
Thuan (1998), subdivided into the metallicity bins covered 
by our SSP model grid, as given in the Zl- and Z2-columns 
in Table 1, for Z = 0.0004 and Z= 0.004 respectively. 

Line ratios in reasonably metal-rich gas 
[Z = 0.008, Zq, and Z= 0.05) are taken from Stasihska 
(1984) as they have been shown to be in good agreement 
with Galactic HII region data (Sivan et al. 1986). No 
further distinction is made between the 3 metallicities 
Z= 0.008, Zq, and Z— 0.05, as galactic HII regions do 
show this full range of metallicities and their line ratios 
at fixed metallicity show considerable scatter. 

The line ratios, of course, do vary with electron tem- 
perature, number density and/or the filling factor of the 
ionized gas. However, we aim to describe typical environ- 
ments at the different metallicities. We have investigated 
the dispersion of line ratios within each of our metallicity 
bins Zl, Z2 and Z3 — Z5 and the impact on integrated 
magnitudes in the broad band filters we consider. We find 
that, on average, the scatter in the line ratios is ^ 30% 
(with few exceptions for weak lines). This scatter trans- 
lates into a scatter in magnitude of at most 0.2 mag for 
solar metallicity, and up to 0.4 mag for the lowest metal- 
licity. Compared with an error of 0.4 mag and 1.0 mag, re- 
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spectively, by not taking gaseous emission into account at 
all, we consider this an improvement. However, we would 
like again to emphasize, that our models are meant to 
improve upon models neglecting the gaseous emission in 
early evolutionary stages in terms of broad band lumi- 
nosities and colors. For star clusters etc. that have de- 
tailed spectral information to determine all their ionising 
parameters, photoionisation codes will allow for a more 
precise individual description of their line strengths and 
their impact on colors. 

In including the gaseous line emission our models go 
beyond Leitherer et al. 's STARBURST99 models which 
do include the gaseous continuum but no emission lines. 
Smith et al. (2002) will soon provide an updated version 
of the STARBURST99 models that also include emission 
lines with particular focus on WR star diagnostics. We 
recall that while STARBURST99 is based on stellar evolu- 
tion input physics from the Geneva group, our models use 
Padova isochrones. 



2.3. The role of white dwarfs 

The inclusion of white dwarfs (WD)s changes the magni- 
tudes of passbands redwards of U by up to few times 0.001 
mag. UV- and FUV-passbands show changes of generally 
0.01 mag. However, even these values are likely to be up- 
per limits - at least in the description of star clusters - 
since their derivation is based on the assumption of the 
presence of a sufficiently large amount of interstellar mat- 
ter (ISM) to be ionized. 

In the case of WDs in star clusters, two aspects are 
of importance: First, WDs only appear after a delay of 
roughly 500 Myr, when SNe had more than enough time 
to blow away the interstellar material from a star clus- 
ter region. Second, the phase of planetary nebulae (PNe), 
when the region surrounding the WD is again refilled 
with gaseous material from the ejected envelopes until the 
ejected material is dispersed into the ISM, is fairly short 
in comparison with the time steps available (~ 25,000 yrs 
vs. 4,000,000 yrs). AS a result, the filling factor (and hence 
the average density) of the ISM during phases when WDs 
could contribute to the emission line flux is expected to 
be well below 1, resulting in a negligible contribution of 
gaseous emission due to WDs to broad band magnitudes. 
In addition, PNe in globular clusters seems to be surpris- 
ingly rare (Jacoby et al. 1997). 

In the case of galaxies, two scenarios might occur: 
Either the galaxy is gas-poor, hence not much matter is 
available that could be ionized, and the situation is similar 
to that of star clusters with low or negligible contributions 
of WDs. Or the galaxy is gas-rich, in which case the emis- 
sion is probably dominated by star formation anyway, as 
described with our models presented here. 



Table 1. Non- hydrogen emission lines and their line 
strengths, normalized to Hp line strength, as a func- 
tion of metallicity (Zl = 0.0004, Z2 = 0.004, Z3 = 0.008, 
Z4 = 0.02 = Z Q , Zh = 0.05). 
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2.4. Filter Systems 

In addition to the Standard Johnson U BVRIJHKL, 
Thuan & Gunn g r, Kron & Koo U+ J+ F+ N+ , 
Washington C M Tl T2, Stromgrcn u v b y, HST WFPC2 
F160BW . . . F814W filters already included in Schulz et 
al. (2002), we here additionally include the HST WFPC2 
medium band filters F410M, F467M and F547M, the 
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NICMOS filters F110W, F160W, F205W and F222M, and 
the ACS/WFC broad band niters F435W . . . F850LP. 

Calibrations for the ACS filters have been performed 
according to Gilliland 2002 (priv. comm.). 



3. Time Evolution of SSP Spectra with Gaseous 
Emission 

The new SSP models presented here extend to younger 
ages and have a better time resolution as compared to 
Schulz et al. (2002). Time resolution now is 4 Myr up to an 
age of 2.35 Gyr, and 20 Myr for larger ages. The youngest 
age available now is 4 Myr as compared to 140 Myr in 
Schulz et al. . In Fig. we present the early evolution 
of the UV - optical parts of SSP spectra of solar and 
low metallicity, respectively, for comparison. The much 
stronger emission lines in the low metallicity spectra are 
clearly seen as well as some differences in line ratios, as e.g. 
that of [Oil] 3727 relative to the adjacent Balmer lines H$ 
and H 7 . The time evolution and metallicity dependence 
of the line strengths will further be discussed in Section 
[J] Fig. |3 shows spectra of a young low metallicity (Z = 
0.0004) SSP over a longer wavelength range. 



4. Impact of Gaseous Emission on Broad Band 
Luminosities and Colors 

4.1. Relative contributions of the gaseous emission - 
lines and continuum - to broad band fluxes 

In Fig. |21 the relative contributions of the gaseous emis- 
sion to broad band fluxes in the Johnson passbands U . 
. . K are shown for some young ages, again for solar and 
low metallicity. The metallicity dependence visible in the 
B-band is due to strong oxygen and neon lines. H^, H 7 
and contribute significantly as well. The V^-band con- 
tribution is dominated by strong oxygen lines (and to a 
lesser extent helium, nitrogen and sulphur), the strengths 
of which are metallicity dependent. The strong contribu- 
tion in the i?-band is caused by H a , which is located close 
to the throughput maximum of the Johnson R filter. The 
/-band is dominated by sulphur lines for high metallici- 
ties (predominantly the [SIII]A9068.60 and [SHI] 9530.85 
lines, both being moderately strong and close to the I- 
band throughput maximum) , which are absent in the case 
of Z = 0.0004. The remaining lines are weak or far from 
the maximum of the filter response function, resulting in 
a prominent dip in Figures I^Jj/Et- Only lines from the 
Paschen- and Brackett-series with small relative flux con- 
tributions are present in the NIR bands J, H and K. 
Here, the contribution from the continuum emission be- 
comes important. 

To interpret the differences between Scalo and Salpeter 
IMF, as seen in Fig. |3 two points are of importance: the 
contributions of the stars to the stellar continuum and 
the contributions to the ionising photons. The stellar con- 
tinuum is dominated by stars around 9 M©, more or less 



500 



Evolution of SSP spectra for Z=0.02 with gaseous emission 
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Fig. 1. Spectrum in terms of flux F\ as a function of wave- 
length A at 3 different times for an SSP of solar metallicity 
(a) and metallicity Z = 0.0004 (b), both with Salpeter 
IMF. 



independent of IMF and metallicity. However, the stars re- 
sponsible for the ionising photons span a larger mass range 
in the low metallicity case as compared to solar metallic- 
ity, providing more ionising photons at low metallicity. In 
combination with the lower overall stellar luminosity for 
the Scalo IMF in comparison to the Salpeter IMF, this 
leads to an enhanced relative emission contribution for 
the low metallicity Scalo IMF, and a reduced contribu- 
tion for the solar metallicity case, both compared to the 
respective Salpeter values. 

While Fig.[3]shows the overall emission contribution of 
the gas to broad band fluxes U . . . K, this is decomposed 
into its line and continuum components in Fig. This 
decomposition of the relative flux contribution of the gas 
to the fluxes in various filters shows that - for all metal- 
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Comparison of 4 Myr spectra with vs. without emission 




wavelength [Angstroem] 



Fig. 2. Spectrum in terms of flux F\ as a function of wave- 
length A for an SSP with Z = 0.0004 and Salpeter IMF 
at an age of 4 Myr, with and without gaseous emission, in 
log-log presentation for better visibility. 



licities - lines dominate in optical filters while continuum 
emission comes into play in the NIR bands. A similar re- 
sult was already reported by Kriiger et al. (1995) in the 
context of evolutionary synthesis models for BCD galax- 
ies. Like the hydrogen lines, the continuum emission is 
not directly metallicity-dependent, but indirectly via the 
higher ionising flux of lower mctallicity SSPs and the dust 
absorption rate f{Z) for ionising photons. 



4.2. Time evolution of gaseous emission contributions 
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Fig. 3. Relative contributions of the gaseous emission - 
lines and continuum - to broad band fluxes U . . . K at 
solar metallicity (a) and low metallicity Z = 0.0004 (b). 



As can be seen from Fig. the gaseous emission con- 
tribution is much stronger as well as longer-lasting at 
lower metallicity. In the Padova (as well as in Geneva) 
isochrones, the most massive stars have longer lifetimes at 
lower metallicities, as compared with solar. As discussed 
in Sect. I2.2I the population of massive stars gets hotter 
and more luminous, on average, at lower metallicities as 
compared to solar resulting in a higher output rate of ion- 
ising photons. In addition, the fraction of ionising photons 
directly absorbed by dust gets lower. Hence the relative 
contribution of gaseous emission e.g. to the i?-band (dom- 
inated by H Q ) starts from 35 % at Zq and from 60 % 
at Z — 0.0004 and decreases by a factor 3 within 5 and 
13 Myr, respectively. At solar metallicity, the overall rel- 
ative flux contribution of the ionised gas falls below 5 % 
in all bands at ages < 8 Myr. At our lowest metallicity, 
Z = 0.0004, it falls below 5 % at 13 Myr (/-band), 15 Myr 
(U, K), and 18 Myr (V, R). 



4.3. Effects of gaseous emission on the color evolution 
of SSPs 

As seen in Fig.H3on the example of V — I the inclusion of 
gaseous emission significantly affects the color evolution 
of SSPs at early stages, depending on the metallicity, and 
the luminosity of a cluster of given mass, as can be seen 
in Fig. [7| This also, of course, affects the age and mass 
determinations from observed colors and luminosities. 

From an observed V —I = —0.2 of a very young cluster, 
e.g., ages of 4 and 7.5-12 Myr would be derived in case of 
solar and metallicity Z — 0.004, respectively, if the effects 
of gaseous emission were not taken into account. Including 
them appropriately, the same observed color V — I = —0.2 
would result in ages 4 and 14.5 Myr, for the two metal- 
licities. This shows that ages estimated on the basis of 
V — I, without taking gaseous emission into account, for 
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Time evolution of gas emission contribution (for Z=0.02) 
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Fig. 4. Decomposition of the total gaseous emission in 
terms of lines and continuum contributions in various fil- 
ters bands U ... K at an age of 4 Myr at solar metallicity 
(a) and low metallicity Z = 0.0004 (b). 
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Fig. 5. Time evolution of the gaseous emission contribu- 
tion to broad band fluxes U, B, V, I, and K at solar 
metallicity (a) and low metallicity Z = 0.0004 (b). 



low metallicity clusters younger than ~ f5 Myr are un- 
derestimated by up to a factor ~ 2. Differences decrease 
with increasing metallicity, and if gas emission is neglected 
at solar metallicity it has become negligible. In a similar 
way, masses estimated from the observed cluster My with 
model mass-to-light ratios that do not account for gaseous 
emission contributions can be overestimated by a factor 2 
at low metallicity. 

Observed colors V— I in the range [-0.35,-0.9] cannot at 
all be accounted for by models neglecting gaseous emission 
contributions. 

5. Discussion 

The aim of our models as presented here is to provide 
the spectral and photometric evolution of SSPs of various 
metallicities that also allow to describe the earliest stages 



of evolution when the ionising flux of a star cluster, a 
stellar (sub-)population etc. is strong enough to cause the 
surrounding gas to emit line and continuum radiation. Our 
motivation was that many HST observations of young star 
cluster systems in interacting and merging galaxies indeed 
do show such very young clusters with broad band colors 
that cannot be interpreted without accounting for their 
gaseous emission. Spectroscopy of some of those bright 
clusters has explicitly revealed their emission lines, e.g. 
Schweizer & Seitzer (1998) in NGC 7252, Colina et al. 
(2002) in NGC 4303, Maoz et al. (2001) in NGC 1512 and 
NGC 5248, ... 

We used most recent compilations of stellar output 
rates in ionising photons. Note that, while finding con- 
siderable differences in the output rates of He-ionising 
photons, the latest and most sophisticated calculations by 
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Fig. 6. Time evolution of the Johnson colors V — I with 
and without gaseous emission at solar metallicity Z — 0.02 
and low metallicity Z = 0.0004. 
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Time evolution of Johnson V-band magnitude with vs. without gas emission 
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Fig. 7. Time evolution of the magnitude in the Johnson 
V-band with and without gaseous emission at solar metal- 
licity Z = 0.02 and low metallicity Z = 0.0004 for SSPs 
of the same mass. 



In Schulz et al. (2002) we gave evolutionary and cos- 
mological corrections for our SSP models without gas for 
convenient use in cosmological structure formation mod- 
els or any kind of models that build up galaxies with ex- 
tended star formation (and possibly also chemical enrich- 
ment) histories. As the changes caused by the inclusion of 
gaseous emission are constrained to the earliest evolution- 
ary stages (t ^ 15 Myr) their effects on the cosmological 
and evolutionary corrections are only visible at redshifts 
extremely close to the redshift of the very onset of star 
formation. We therefore do not present cosmological and 
evolutionary corrections here and advise potential users 
to stay with the tables in Schulz et al. (2002). 

Comparison of our GALEV model using Padova stellar 
input physics with results from Leitherer et al. 's (1999) 
STARBURST99 code, using Geneva stellar tracks, during 
the first Gyr of evolution covered by the latter shows fairly 
good agreement. As an example we show in Fig. [S] the 
time evolution of B — V . The differences are minor (< 
0.1 mag), and mainly due to the better time resolution of 
the STARBURST99 code. Comparison of colors including 
longer wavelengths, e.g. V — I or V — K, show larger de- 
viations that are due to different input physics, especially 
a different treatment of red supergiants and the omission 
of the TP-AGB phase in starburst99. 
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0.25 



Comparison of our code with STARBURST99, solar metallicity 




GALEV, with EL 
GALEV, without EL 
SB99, standard mass-loss 
SB99, high mass-loss 
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Smith et al. (2002) confirm the rates from Schaerer & de 
Koter (1997) that we use for H-ionising photons. 

Our focus is on the impact of gaseous emission on the 
early luminosity and color evolution. For low metallicities 
(Z < 0.004), we chose to use line ratios for the different 
metallicities of our SSP models as observed in HII regions 
of the respective metallicities. The line ratios for higher 
metallicities were taken from the theoretical models pre- 
sented by Stasihska (1984). Hence, our models should not 
be used for any analyses of line ratios. For this purpose 
we recommend the use of models that couple to photoion- 
ization codes like those of Chariot & Longhetti (2001) or 
Moy et al. (2001). 



Fig. 8. Comparison of our models with STARBURST99 
models with different mass-loss rates, for solar metallic- 
ity and the Johnson color B — V. 



6. Structure and Description of Electronic 
Data-files 



Our model results are accessible at 



http://www.uni-sw.gwdg.de/~galev/panders/ 



CDS and at 
The files 



have the following naming conventions: 



8 



P. Anders & U. Fritze - v. Alvensleben: Young Stellar Populations 



— Each directory is labelled ssp_model_zx_y. Here zx (x 
= 1 ... 5) indicates the metallicity (zl <^ Z = 0.0004, 
z2 & Z = 0.004, z3 & Z = 0.008, z4 <S> Z = 0.02, 
z5 <^> Z = 0.05), and y the assumed IMF $(m) ~ 
m Q ('salpeter' <^> a = —2.35 for all masses; 'scalo' <^> 
a = -1.25 for m < 1M , a = -2.35 for 1M© < m < 
2M Q and a = -3.00 for 2M Q < m). 

— Each directory contains a file named 'plotsgal.dat' con- 
taining the cluster mass as a function of age, a file 
named 'galspec.dat' containing the integrated cluster 
spectra as a function of age, and a file named 'mag- 
nitudes. dat' containing the integrated magnitudes of a 
cluster as a function of age in various passbands. 

— Further information is included in the README file. 

7. Conclusions and Outlook 

We present an update of the evolutionary spectral syn- 
thesis models for SSPs of various metallicities presented 
in Schulz et al. (2002). Using the same stellar isochrones 
provided by the Padova group, including the TP-AGB 
phase important for V — I and V — K colors of SSPs 
in the age range 10 8 . . . 10 9 yr, we add the gaseous 
emission contributions, both lines and continuum. They 
are very important not only for the spectral evolution 
but also for the photometric evolution in terms of broad 
band colors in early evolutionary stages, as long as ion- 
ising stars are alive. We use the best available output 
rates for H-ionising photons and observationally sup- 
ported emission line ratios appropriate for the different 
metallicities. This allows us to extend our SSP models to- 
wards younger ages than accessible before. We provide 
an extensive set of electronic data for the time evolu- 
tion of UV through NIR spectra of SSPs of metallic- 
ities Z = 0.0004, 0.004, 0.008, 0.02 =Z Q , and 0.05, for 
Salpeter and Scalo IMFs, from ages as young as 4 Myr 
all through 14 Gyr in timesteps of 4 Myr until an age 
of 2.35 Gyr and of 20 Myr thereafter. We also provide 
the luminosity and color evolutions in a large number of 
ground- and space-based filter systems including Johnson, 
Washington, HST WFPC2, NICMOS and ACS/WFC. 

Our results also allow for easy superposition of SSPs 
of various ages - and possibly metallicities - to describe 
galaxies with all possible star formation and metal en- 
richment histories. They also are readily included into dy- 
namical and cosmological galaxy or structure formation 
scenarios that contain a star formation criterium. 

The evolutionary and cosmological corrections we gave 
in Schulz et al. (2002) are not affected by gaseous emis- 
sion as it is constrained to stages earlier than 15 Myr, i.e. 
to redshifts extremely close to the redshift at which star 
formation sets in. 
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